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Abstract: Using the twisted partition function on M'^ x S^, we argue that the deconfine- 
ment phase transition in pure Yang-Mills theory for all simple gauge groups is continuously 
connected to a quantum phase transition that can be studied in a controlled way. We explic- 
itly consider two classes of theories, gauge theories with a center symmetry, such as SU{Nc) 
gauge theory for arbitrary Nc, and theories without a center symmetry, such as G2 gauge 
theory. The mechanism governing the phase transition is universal and valid for all simple 
groups. The perturbative one-loop potential as well as monopole-instantons generate attrac- 
tion among the eigenvalues of the Wilson line. This is counter-acted by neutral bions — 
topological excitations which generate eigenvalue repulsion for all simple groups. The transi- 
tion is driven by the competition between these three effects. We study the transition in more 
detail for the gauge groups SU{Nc), Nc > 3, and G2. In the case of G2 there is no change 
of symmetry, but the expectation value of the Wilson line exhibits a discontinuity. We also 
examine the effect of the 0-angle on the phase transition and critical temperature Tc{6). The 
critical temperature is a multi-branched function, which has a minimum at ^ = vr as a result 
of topological interference. 
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1. Introduction 

In our recent work, we studied a semi-classical mechanism for the center symmetry changing 
phase transition in four dimensional non-abelian SU{2) gauge theory [1]. The main idea of 
that work was to use continuity and the notion of a twisted partition function. We argued 
that the deconfining phase transition of pure Yang-Mills theory on x which takes 
place at strong coupling, is continuously connected to a quantum phase transition that can 
be studied reliably, albeit using non-perturbative methods, at weak coupling. 

The most important non-perturbative effect that arises in the weak coupling realization 
of the center symmetry changing phase transition is governed by a new class of topological 
objects, called "neutral bions". These objects are related to four dimensional instantons, but 
their physical effects are quite different. The 't Hooft vertex for a neutral bion induces a 
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Table 1: The simple Lie algebras g together with their associated compact simply-connected 
Lie groups G and the compact adjoint Lie groups Gad- The last column lists the center 
symmetry of G, which is isomorphic to the fundamental group of Gad- 

center-stabilizing potential in SU{2) theory, by generating a repulsion among the eigenvalues 
of the Wilson line.^ This effect counter-acts the center-destabilizing perturbative potential 
for the Wilson line, which leads to an attraction among the eigenvalues [3]. In our previous 
work we showed that the monopole-instanton amplitudes also lead to an attraction among the 
eigenvalues of Wilson line. We systematically studied the competition between these three 
effects in the case of SU{2) gauge theory. We observed that the twisted partition function on 
X §^ exhibits a second order phase transition, consistent with the second order transition 
observed in lattice gauge simulations of thermal SU{2) Yang Mills theory. 
In the present work, our goal is two-fold: 

aj to generalize the phase transition mechanism in SU{2) to arbitrary compact simply- 
connected Lie groups, G, including the theories without a center, G2,Fa,E8- 

b) to understand the topological 6'-angle dependence of the critical temperature Tc. 

For convenience, we list the simple Lie algebras q together with their associated compact 
simply-connected Lie groups G and the compact adjoint Lie groups Gad iii Table 1. The last 
column lists the center group of G, which is isomorphic to the fundamental group of Gad- For 
theories with adjoint matter, the center symmetry is the same as the center group. 

1.1 Twisted partition function and continuity of the deconfinement transition 
between weak and strong coupling 

We consider Yang-Mills (YM) theory with one adjoint Weyl fermion with mass m. In the 
chiral limit this theory reduces to M = 1 supersymmetric YM theory, and in the limit m — )• 00 

'^For this reason, we sometimes refer to the neutral bion in theories with a center symmetry as "center- 
stabilizing bion" [2]. 
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Figure 1 : Phase diagram of Yang Mills (YM) theory with an arbitrary simply-connected gauge group 
G and a single Weyl fermion of mass m in the adjoint representation, subject to periodic (left) and 
anti-periodic (right) boundary conditions on M.^ x S^. Right: For the ordinary partition function the 
transition occurs at a critical temperature of order A, and is not accessible by analytical means. Left: 
In the case of the twisted partition function, the thermal deconfinement phase transition in pure Yang 
Mills theory, corresponding to the limit m — > oo, is conjectured to be continuously connected to a 
quantum phase transition in supersymmetric Yang Mills (SYM) theory deformed by a gluino mass 
term. For small m, the transition takes place as small L and is analytically tractable. For theories 
without center symmetry, there is no symmetry breaking associated with the transition. 



it corresponds to pure YM theory. Let H and F denote the Hamiltonian and fermion number 
operator, respectively, and L the circumference of the circle. We introduce the twisted 
(non-thermal) partition function: 

Z{L, m) = tr [e-^^(-l)^] = tr [g-^^+^-^j , (LI) 

which, when viewed as a function of fermion mass, interpolates between the supersymmetric 
(Witten) index and the ordinary thermal partition function of the pure Yang-Mills theory. 
Namely, 

Z{L,m = 0) = Is{G), 

Z{L, m = oo) = Z™{^), L = (3 , (L2) 

where IsiG) = hy is the index, which is equal to an invariant /i^ (independent of L) associated 
with the group G and called the dual Coxeter number {h^ equals N^, for SU{Nc)). Since 
dZ{L,0) /dL = for any L, there is no phase transition at m = [4-6]. Z^^(/3) is the 
partition function of the pure Yang-Mills theory. Eqs. (LI) and (1.2) permit us to continue the 
deconfinement phase transition in pure YM theory to a semi-classically calculable transition 
in the small L-m regime as shown in Fig. L 
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1.2 Review of the dynamics for m = and small-L 

The dynamics of A/" = 1 SYM theory with gauge group G in the small-L weak coupling regime 
is that of the abelianized theory, U{lY , where r = rank(g) is the rank of the Lie algebra of 
G [4-6]. The abelian description is valid at distance scales > Lr/27r [7]. This follows from 
the fact that the Wilson line holonomy behaves as a compact "adjoint Higgs field" in the 
weak coupling regime. We note that there is no such description at strong coupling, see the 
detailed discussion in Section 5. In a Euclidean context, the long distance theory is described 
as a dilute gas of topological 1-defects (item (i) below) and 2-defects (items (ii) and (iii)): 

(i) monopole-instantons Mi with i = 1, . . . , rank(g) + 1, 

(ii) magnetic bions Bij = [TWiA^j], V Aij < 0, 

(iii) neutral bions Bu = [MiMi], V An > 0. 

Here, Ay is the extended Cartan matrix of G. Small 4d instantons are responsible for breaking 
the anomalous axial U{1)a symmetry down to an anomaly free discrete subgroup, TL^hy ■ But 
apart from that, they have no sizable impact to any physical phenomena in this regime. In the 
small-L regime, 4d instantons can be viewed as composites of the monopole-instantons, i.e., 
T^Ad ~ 11^=0 [-^i]'^'' ' ■^l^sre A:J are dual Kac-labels (or co-marks; equal to unity for SU{Nc)) 
and Yl^j=o ~ ~ ^s{G) is the dual Coxeter number. This formula embodies the physics 
of fractionalization of a large 4d instanton > Lr into kj monopole-instantons of type Mj, 
each of which carry two adjoint fermionic zero-modes.^ 

Since large 4d-instantons do not exist in this regime the problem of the breakdown of 
the dilute instanton gas approximation, which follows from the fact that instanton size is a 
modulus, is avoided. Moreover, since the 1-defects and 2-defects obey a separation of scales, 

rm <C rb < dm-m < c^b-b, (1-3) 

the treatment of monopole-instantons (Mi) and topological molecules (Bij) in the dilute- 
gas approximation and the effective theory derived from their proliferation in the vacuum is 
reliable. This is depicted in Fig. 2: and rb are the typical sizes of a monopole-instanton 
and bion events, and dm-m and db-b are the typical inter-monopole-instanton and inter-bion 
separation. For a detailed discussion of the scales, see [14]. 

Monopole-instantons have two fermionic zero modes associated with a vertex of the form 
of , ipii^j -H- ^fcA^fc, where W{^i) is the superpotential of the effective theory. This 
means that the superpotential can be extracted from the monopole amplitude, but it does 
not imply that the mass gap for bosonic fluctuations is generated by monopoles, as is evident 
from the existence of a fermionic bilinear in the amplitude. The bosonic potential can be 
computed from the superpotential, V = if^"^ ^ ^ij^ij^ but the physical mechanism 

^The topological classification of finite action solutions on xS'^ can be found in [3]. Monopole constituents 
of the instanton (caloron) solution for non-trivial holonomy were found in [8,9], and the relevant index theorem 
is discussed in [10,11]. The monopole vertex for general gauge group can be found in [6]. Magnetic bions are 
studied in [12,13], and neutral bions are discussed in [1,2,14]. 
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Figure 2: The Euclidean vacuum of the small-m, L theory can be described as a plasma of monopole- 
instantons (gray circles) and anti-monopole-instantons (white circles) with fermionic zero modes (un- 
paired arrows). The paired events are magnetic and neutral bions. Neutral bion amplitudes generate 
repulsion among the eigenvalues of the Wilson line and magnetic bions generate a mass gap for gauge 
fluctuations via a generalization of the Polyakov mechanism to a locally 4d theory. 



which generates the bosonic potential and the mass gap for bosonic fluctuations is due to 
bions, correlated monopole-anti-monopole instantons without any fermionic zero modes. The 
bosonic potential also has /i^ minima, leading, at weak coupling, to the spontaneous breaking 
of the discrete chiral symmetry, TLihy ^2- This, in turn, generates a dynamical mass for 
fermions. The importance of this point of view, apart from providing the correct interpretation 
of the physical phenomena governing the dynamics in the supersymmetric theory, is that semi- 
classical monopole and bion amplitudes also exist in non-supersymmetric theories, where the 
bosonic potential cannot be extracted from the super-potential [2,12-14]. 

1.3 Phase transition in the small m-L regime and universal aspects 

There are two main effects of adding a small fermion mass term. The mass term lifts the zero 
modes of the monopole-instantons. This implies that there is a non-zero monopole-instanton 
contribution to the bosonic potential. The mass term also breaks supersymmetry, which leads 
to a perturbative contribution to the potential for the holonomy [3]. Studying the competition 
between these effects and the bion induced potential already present at m = shows that 
there is a phase transition at some critical compactification scale that grows with m. We find 
a description of this phase transition valid for all Lie groups, G: 

1. Neutral bions always generate repulsion among the eigenvalues of the Wilson line around 
S^. For theories with a Zjv center symmetry, the repulsion leads to a Z jv-symmetric 
distribution, while for theories without a center symmetry, it leads to a non-degenerate 
distribution of eigenvalues, as we show explicitly for G2. 
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2. Monopole-instantons, along with the perturbative potential for the Wilson line, lead to 
an attractive interaction between the eigenvalues of the Wilson line. This interaction 
increases with m. Monopoles prefer configurations in which the eigenvalues accumulate. 
This leads to center-instability whenever Z{G) is non-trivial and to eigenvalue bunching 
whenever Z{G) is trivial. 

3. Magnetic bions lead to an attraction among the eigenvalues of the Wilson line. They also 
generate a mass gap for gauge fluctuations, and they are responsible for the confinement 
of electric charges. However, the combined effect of neutral and magnetic bions (which 
are of the same order in the semi-classical expansion) always generates a repulsion 
among eigenvalues. 

4. The competition of neutral and magnetic bions on the one hand and monopole-instanton 
and perturbative effects on the other hand is responsible for the semi-classical realization 
of deconfinement. For all gauge groups but G2, -F4, -Eg, this phase transition is associated 
with a change in the center symmetry realization. For SU{Nc) with Nc > 3, we find a 
first order phase transition associated with a change in the center symmetry realization. 
In earlier work, we found a second order phase transition for Nc = 2 [1]. In theories 
without a center symmetry, we find a first order transition accompanied by a jump in 
the Polyakov loop. These findings agree with lattice results (see [15] and [16] for a 
recent review), providing support for the conjectured continuity of the deconfinement 
transition depicted in Fig. 1. 

Finally, we investigate the topological 6'-angle dependence of the twisted partition func- 
tion. As is well-known, turning on a 0-angle introduces a sign problem in the Euclidean 
path integral formulation. The semi- classical manifestation of the sign problem is a complex 
fugacity of monopole-instantons. The phase of the amplitude generates interference between 
Euclidean path histories, which can be interpreted as the analytic continuation of Aharonov- 
Bohm type interference. We refer to this phenomenon as "topological interference". In the 
case of finite rank gauge groups topological interference leads to 6 dependence of the critical 
temperature. The minimum Tc occurs at = vr. At this point the 9 dependence is non- 
analytic, and all theories have two degenerate vacua in the confined phase. At = 00 the 
critical temperature is independent of 6. 



2. SU(Nc), Nc > 3: First order center symmetry changing phase transition 

As reviewed in Section 1.2, at sufficiently small L the supersymmetric SU{Nc) theory on 
M.^ X dynamically abelianizes down to C/(l)^'=~^ at distances larger than LNc/2t:. In this 
regime the light bosonic fields are the A'^c — 1 dual photon fields a {a Nc — 1-dimensional 
vector) and their scalar superpartners (p, which are the "uneaten" components of the gauge 
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connection along S^. It is convenient to expand the fields (p and a as follows: 

^ 96 27rk + e ^ ^, , , 

The primed fields describe the fiuctuations of the fields (p and a around their supersymmetric 
and center-symmetric ground state. For every Nc the supersymmetric ground state is labeled 
by an integer k = l,...Nc, the vacuum angle 9, and the Weyl vector p. The integer k 
corresponds to a choice of vacuum corresponding to the spontaneous breaking of the anomaly- 
free discrete chiral {R-) symmetry, Z2Ar^ Z2. 

We normalize the SU{Nc) Cartan generators H = {Hi, . . . ^Hn^^i) by tr:[HaHh] = Sab- 
All the SU{Nc) roots have length squared equal to 2, and the roots and co-roots, a and a*, 
as well as the fundamental weights and co- weights, Co and oj* , are the same. The Weyl vector 
is given by p = The fundamental weights obey Wj • dj = 6ij = 1, . . . Nc — I), and 

the Weyl vector obeys p- Si = 1, foic i = 1, . . . Nc — 1, and p- Sn^ = 1 — Nc- 

The fields b' and a' in (2.1) are periodic, 

b' r^b' + %-uj, a r^a + Ittuj . (2.2) 

Note that for small g periodicity of b' is irrelevant. The gauge holonomy (Wilson loop around 
S^) in terms of the fields (2.1) is given by: 

J7 = exp A —H ■ p + i —H ■ b'\ . (2.3) 

At 6' = 0, obeys tr Vt^ = 0, V n 7^ (mod Nc). We note that tr Jl" are order parameters of 
the Zjv^ global center symmetry that the adjoint-fermion SU{Nc) theory acquires when com- 
pactified on x and that center symmetry is unbroken in any of the Nc supersymmetric 
vacua. 

It is also important to note that the Wilson line that satisfies (tr 0") = at weak 
coupling is associated with the eigenvalue configuration depicted on Figure 3(b). We refer 
to this type of gauge holonomy as weak coupling non-trivial holonomy. At strong coupling, 
there is no adjoint Higgsing, as the fiuctuations of the eigenvalues are larger than the typical 
inter-eigenvalue separation, hence in the latter, the eigenvalues are essentially randomized 
over the dual circle on which the eigenvalues live. We will discuss this distinction in more 
detail in Section 5. 

2.1 Relevant scales and two large-Nc limits: 't Hooft vs. Abelian 

The gauge theory compactified on x contains an extra parameter, the compactification 
scale L. Thus, if center symmetry remains unbroken at small L, the large- A'^c limit is more 



-7- 





(a) (b) (c) 



Figure 3: Three types of gauge holonomy: a): Weak coupling trivial holonomy associated with 
deconfined phase, b): Weak coupling non-trivial holonomy associated with the weak coupling confined 
phase, c): Strong coupling non-trivial holonomy associated with the strong coupling confined phase. 
It is important to note that b) admits an adjoint-Higgs description, whereas c) does not. These two 
confinement regimes are continuously connected. 

subtle than the standard 't Hooft limit. ^ In fact, there are two large- A'^c limits, which we refer 
to as the (naive) 't Hooft limit and the abehan limit. The reason that there is more than 
one limit is that mw, the mass of the lightest off-diagonal gluon, goes to zero if the large- A'^c 
limit is taken with LA fixed. This means that the effective abelian theory breaks down in the 
naive 't Hooft hmit. When we study the large- A'^c limit we will therefore consider the limit 
Nc ^ oo with niw = 2it/{NcL) = const, i.e., L = 0{N^^). In other words, ensuring the 
validity of the weakly-coupled effective theory requires shrinking the size of the circle when 
taking the large- A^c limit. In what follows, we will express all dimensionful quantities in the 
effective lagrangian for the light modes on x in terms of and the scale parameter 
of the gauge theory. A, with the understanding that they are both fixed in the large- A'c limit, 

^ = /U exp 



m^ = —»A. (2.4) 

The above expression for A is the standard 2-loop expression for the scale parameter [17]. 
Note that Ref. [6] uses a definition of A'^ that differs by a factor of proportional to Nc- The 
gaugino mass m introduces another length scale in the problem. The validity of the fermion 
zero mode induced pairing mechanism of bions requires that m <^ g'^{myi/) / L, the inverse 
bion size. In terms of the scales in Eq. (2.4), the gaugino mass should parametrically obey 

- « ^ ^ (2-5) 

^This issue does not arise in thermal compactifications, and as a consequence it has not been discussed 
until recently, see [7]. 
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an upper bound which remains finite in the large- A'^c hmit. 



2.2 The bion and monopole potentials in softly-broken SU(Nc) theory on X §1 

Here, we describe the Lagrangian governing the long-distance dynamics of the SU{Nc) gauge 
theories on x for general Nc- The effective Lagrangian follows from the superpotential 
of Ref. [6] and the field redefinition in Eq. (2.1) in a rather straightforward way, generalizing 
the detailed derivation given in [1] for Nc = 2. The kinetic term of the light fields (2.1) is:^ 



kin. 



167r2L 



(2.6) 



where mw is given in (2.4). The nonperturbative potential due to bions and monopole- 
instantons is: 



np 



^bion 



Nc 



i=l 
- Nc 

E 



-2d,-b' 



[( 



1 + 



Si ■ b' 



e cos 



a,; - Oj+i) • (?'] 

2TTk + 9 



(2.7) 



a,; • fJ + 



where we have included the dependence on the 9 parameter of the gauge group, which only 
contributes if the gaugino mass is nonzero. The parameter k = 0,l,...,A'^c — 1 labels the 
vacuum branch as in Eq. (2.1). The true effective potential corresponds to the branch that 
gives the minimum ground state energy, the other branches describe metastable vacua. These 
vacua are expected to become quasi-stable in the large- A'^c limit [18]. The strengths of the 
bion and monopole potentials are given by 





bion 



exp 



167r2 



exp 



87r2 
9'(/")iVc 



(2.8) 
(2.9) 



IKf^gKf^Nc 

For Nc = 2 and 9 = the above result reduces to the expressions given in [1]. We can express 
these results in terms of the physical scales defined in Eq. (2.4). The kinetic term is 

^ 1 mw 
127rlog(^ 

and the coefficients (2.8, 2.9) of the bion and monopole potentials are 

27 . .. ^ 



(2.10) 



^bion 



9 mA^ 

2-K my/ 



(mw\ 
log ) , 



log . 



V A y 



(2.11) 
(2.12) 



''The kinetic terms of a and a' receive perturbative corrections along the Coulomb branch. These are closely 
related to the non-cancelling one-loop fluctuation determinants around supersymmetric instanton-monopole 
solutions on 'M? x S'. This subtlety is studied in detail in Appendix A of [1]. However, as in that reference, 
our results in the weakly-coupled calculable regime are unaffected to leading order. 
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For future use, we also introduce the dimensionless ratio of Eq. (2.12) and (2.11): 



yLn~ 3A3 -3A(ALiV,)2- ^'-'^^ 

We first study the spectrum in the hmit when the supersymmetry breaking perturbation is 
turned off, m = 0. For this purpose we expand the bion potential to quadratic order in b' 
and a' . Because of the exact supersymmetry the b' and a' masses are identical and we can 
concentrate on the b' field. We rescale b' to canonically normalize the kinetic term (2.10), 
and also drop the prime in the following. We find the following quadratic lagrangian 



1 

£ = -{d,bf + mlY^ 



i=l 

with 



Si-b] - (di-b) ( cti+i ■ b 



(2.14) 



The mass term in (2.14) can be diagonalized by switching to U{N^ roots given by N^- 
dimensional vectors of the form ai = (1,-1,0,...), a2 = (0,1,-1,0,...), etc., oat^ = 
(— 1, 0, . . . , 0, 1), see [19]. This introduces an extra massless particle that decouples from 
the rest of the spectrum. We get 

1 1 ^= 

^ = l^idi^bf + - mlY,{hi+2 - 2bi+i + bif . (2.16) 

i=l 

This is a simple lattice model that can be diagonalized by introducing Z^v^ Fourier modes 



p=0 



The spectrum of the b fields is 

m, = 4mosin2(^)=18A^^^^^sin2(^) , (p = 1, . . . , iV, - 1) , (2.18) 

where we have dropped the unphysical massless mode p = (which was introduced above via 
the transition to A'^c- dimensional roots). From Eq. (2.18) we conclude that all the 6-fields are 
parametrically lighter than the VF-bosons and that the effective theory is consistent in the 
large- A^c hmit, see Eq. (2.4). 

Let us now turn on the gaugino mass m. At first, we ignore the ^-dependence and take 
= and A; = (it is easy to see that for = the /c = vacuum in Eq. (2.1) has the minimal 
energy). In this vacuum, there is no mixing between the b' and a' fields when the monopole 
potential is expanded up to quadratic terms around their vanishing expectation values. In 
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the weak coupling limit = [3 log(mvy/A)] ^, so we can neglect the g'^Nc term in the 

monopole potential in Eq. (2.7). We now find, instead of (2.14), the quadratic Lagrangian 



C 



i=l 



di ■ b 



ai-b] ai+i ■ b 



-{d^■bf 



(2.19) 



Thus, the monopole term is diagonal in the Fourier basis given in Eq. (2.17), and the only 
difference is that the quadratic part of the monopole potential has eigenvalues that are pro- 
portional to the square of sin(p7r/A^c)- The squared mass of the b field as a function of m 
is 



ml 



sm 



Cm 

T 



sm 



(2.20) 



We observe that the center symmetric vacuum becomes locally unstable for Cm > c^, where. 



taking p 



1 c* 



4sin (-^). This corresponds to m > m* with 



3A3 



m 



sm 



m 



W 



(2.21) 



for fixed L as one varies m. On the other hand, for fixed m, as one varies L, this implies a 
local instability at 

/ Am 



L* 



A- 



1 



Nc 



A' 



4m 
3A 



(2.22) 



It is clear that m* (or L* ) is within the region of validity of the semi-classical approximation 
given in Eq. (2.5). 

We note, however, that for any > 2 the non-perturbative potential contains cubic 
terms. This implies that the transition is first order, and that is not the critical coupling, 
but rather it is the limit of metastability above which the confining vacuum ceases to be a 
(local or a global) minimum. The critical mass, rn-cr, is smaller than m* in Eq. (2.21). We will 
determine the critical values of in Section 2.4. For < Cm < c*^ the center-symmetric 
minimum is a metastable local minimum but not the global minimum. 



2.3 The perturbative potential for arbitrary Nc 

Before we continue our study of the phase transition we will first show that the perturbative 
(Gross-Pisarski-Yaffe) contribution to the potential for the holonomy is sub-leading compared 
to the monopole and bion induced potentials, as in the Nc = 2 case. The perturbative 
potential for the Polyakov line in SU (Nc) is 

,2 



Vr 



pert 



m 



Nc>j>i>l 



Oil 

2 



(2.23) 



where we have neglected higher order terms in the mass of the adjoint fermion. Here, B2{x) 



X 



x + 1/6 is the Bernoulli polynomial of order 2, and 



(2.24) 



- 11 - 



For small b[ we can ignore the periodicity of the potential. The center symmetric point 6^ = 
is a local maximum of the perturbative one-loop potential. At order m? the potential near 
the center symmetric vacuum is exactly quadratic, 



m 



pert 



pert 

2N7 



m 



j>i 



TVW 

A 



pert 31og(21Wl) 



(2.25) 



where we have dropped a constant term and normalized the kinetic term as in Eq. (2.14). 
The sum can be diagonalized in terms of the Z^r^ Fourier modes given in Eq. (2.17). The 
tachyonic mass for the p'th mode is 
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independent of p. We can compare this result to bion generated mass of the lowest mode, 
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At the critical mass m* defined in equ. (2.21) the tachyonic mass is 
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(2.28) 



and we conclude that the perturbative contribution is suppressed by three powers of \og{my/ / A) . 
This is the same suppression factor we found in A^c = 2 [1]. We will henceforth ignore the 
contribution of the perturbative fluctuations to the potential for the holonomy. 

2.4 The phase transition and the Hmits of metastabiUty for Nc > 3 

In Section 2.2 we determined the limit of metastability at which the origin of the effective 
potential for Nc> 3 becomes unstable, and the metastable phase disappears. In this section 
we will numerically determine the critical for the first order phase transition. 

SU(3): For Nc = 3, one can plot the potential as a function of the deviation of the holonomies 
from their center symmetric values (i.e., &X'^2) ^'-'^ different values of Cm and observe the 
nature of the transition, see Figure 4. For 6 = 0, the expectation value of a' remains zero. 
The critical value for the first order transition is = 2.446, see Figure 5. 

For Cm < c^, where is one of the limits of metastability, there is a unique center- 
symmetric minimum (confining phase), see the left panel of Figure 4 as well as Figure 5. 
For Cm < Cm < Cm, there is a center-symmetric global minimum (confined phase) and three 
center-broken metastable (deconfined) local minima (not shown in Figure 4, but shown in 
Figure 5). This also provides an alternative definition c!^. It is the value of Cm below which 
the metastable center-broken minima disappears, hence the name limit of meta-stability. 

For c'f^ < Cm < Cm, where cj^ = 3 is the limit of metastability, the three Z3 breaking 
minima are global degenerate minima, and the center-symmetric confining phase is a local but 
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-1.0 -0.5 0.0 0.5 1.0 -1.0 -0.5 0.0 0.5 1.0 

Figure 4: Contour plots of the bion- and monopole-instanton-induced potential, as a function of the 
two holonomies, showing the first order phase transition for SU(3) (darker shades represent smaller 
values of the potential). Left panel: Contour plot for Cm < cj^ < (cm = 2.20, = 2.446) as 
a function of 61,62- The Zs-symmetric (confining) minimum is at the origin. Right panel: Contour 
plot for < c,„ < c*j (cm = 2.5) as a function of 61, &2- The Zs-breaking global minima are clearly 
visible, and the Za-symmetric confining minimum is meta-stable. 

not a global minimum, see the right panel of Figure 4. In this regime the confining phase is 
meta-stable. Finally, for Cm > c^, the center-symmetric point ceases to be a local minimum, 
and this correspond to the other limit of metastability. This case is not shown in Figure 4, 
but shown in Figure 5. 

SU(Nc),Nc > 3: The general structure that emerges for SU{Nc),Nc > 3 is similar to 
the SU{3) case shown in Figure 5. We have four characteristic domains for the bion and 
monopole-instanton induced potential: 

• Cm < or L > L**: There is a unique center-symmetric (confined) minimum. 

• < Cm < or L** > L > L^'^: A global center-symmetric (confined) minimum and 
Nc meta-stable Z^r^ breaking (deconfined) minima. 

• Cm < Cm. < Cm or L^^' > L > L* : A metastable center-symmetric (confined) minimum 
and Nc global Z^v^ breaking (deconfined) minima. 

• Cm < Cm or L* > L: Nc center-breaking global (deconfined) minima. 

A lattice study of the endpoint of the regime of metastability in pure Yang-Mills theory 
was reported in [20]. This study was motivated by the old idea that in large- A'^c QCD the 
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Figure 5: Effective potential due to bions and monopoles for different values of Cm- The potential 
(in units of V^^^^) is shown as a function of bi for 62 = 0, corresponding to a cut along the a;-axis in 
Fig. 4. 

Hagedorn temperature Tjj, the limiting temperature of a hadronic resonance gas, must cor- 
respond to a second order phase transition which occurs above the (first order) deconfining 
transition [21-23]. The regime Tc < T < Th is therefore metastable. This behavior was also 
observed in large- iVc theories compactified on on x 5^ [24,25]. 

For A'^c > 3 the potential is a function of A'^c — 1 holonomies. Instead of trying to plot 
this function we will study the vacuum expectation values of the Nc eigenvalues of the Wilson 
line around as a function of the gaugino mass parameter Cm- We will, again, observe a 
discontinuous transition, which occurs at values of smaller than the cj^ given in (2.21). 
The eigenvalue repulsion due to the neutral (center-stabilizing) bions is dominant at small 
Cm and is, upon increasing Cm, countered by the eigenvalue attraction due to monopole- 
instantons along with the one-loop potential which is sub-leading, see Section 2.3. We plot 
the eigenvalues of Q. crossed. We show the behavior of the eigenvalue distributions for 

Nc = 4, 5, and 10 in Figures 6, 7, and 8, respectively. In order to make these plots we chose 
a specific value of the 't Hooft coupling, g^Nc = 0.1. The critical value is independent 
of the coupling constant in the weak coupling limit. Results for 0*^^^ and are compiled in 
Table 2. 

The discontinuity at the first order phase transition can be quantified in terms of the 
change in the trace of the Wilson line. In weak coupling, g'^Nc/{4:7r) <C 1, we have 

tr{n) ^ i (^^) i-tr {n,H ■ {h')) , (2.29) 

where VLq = e^^'^ is the center-symmetric holonomy. The discontinuity A|tr(r2)| at c^, 
calculated from Eq. (2.29) using the numerical value for c^, is also given in Table 2. The 
data for the discontinuity of the trace (with maximal trace normalized to unity) can be 
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1.382 


0.97 


0.61 


0.75 


6 


1.000 


0.68 


0.62 


0.91 


7 


0.753 


0.50 


0.62 


1.07 


8 


0.586 


0.38 


0.61 


1.22 


9 


0.467 


0.30 


0.61 


1.37 


10 


0.382 


0.24 


0.61 


1.53 



Table 2: Critical values cj^ and for the endpoint of the metastable phase and the first 
order transition in SU{Nc) gauge theory with one adjoint fermion on x Si. We also show 
the critical coupling scaled by Att'^/N^, and give the discontinuity of the Wilson line in 
the weak coupling limit. 



(cl) Cm (c) Cm — 

Figure 6: The distribution of the eigenvalues of the Polyakov loop around S}^ for Nc = A for different 
values of Cm, shown for g^A^c = 0.1. 



represented by the expression 



trn 



0.163 



0.085 g'^Nc 



(2.30) 



which is valid in the large- A'^c limit with the parameters given in Eq. (2.4) held fixed. 



2.5 More on the (abelian) large-Nc limit and (evading) the Hagedorn instability 

So far, we have computed the critical mass numerically for Nc = (3,4, . . .). The numerical 
results indicate that the critical value of the gaugino mass, nicr, is equal to up to a factor 
of order one for arbitrary Nc- This can be seen from the fourth column of Table 2, which 
shows the value of scaled by the asymptotic value of a^. Therefore, using Eq. (2.21, 2.22), 



we observe that 



^ ~ liLAf at fixed L, or L* ~ A"Ym ^^^"^ conclude that in 

the scaling regime the SU{2) result derived in [1] has a smooth large- A^c limit. 
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(a) Cm ^ c„ 



(b) c„ 



0.97 



Figure 7: The distribution of the eigenvalues of the Polyakov loop around Ei\ for Nc = 5 for different 
values of c„j, shown for g^N^ — 0.1. 
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Figure 8: The distribution of the eigenvalues of the Polyakov loop around E>]^ for Nc — 10 for different 
values of c^, shown for gr^iVc = 0.1. 



We note that in the abelian confinement regime the size of the dual circle S on which the 
eigenvalues reside is equal to 1/L, which grows linearly with Nc, while the separation of the 
eigenvalues remains fixed, 2n/{LNc) ~ 0{N^). This is in sharp contrast with the non-abelian 
confinement regime and the ordinary large- A'^c limit. In the latter, 1/L = 0{N^) and the 
separation between eigenvalues is 2tt/{LNc) ~ 0{N~^), forming a dense set in perturbation 
theory. In the latter case, since 2tt/{LNc) <C A, all the low momentum modes are strongly 
coupled, and the eigenvalues are uniformly distributed over the unit circle. At the critical 
point in the weakly coupled regime, the uniform separation between eigenvalues exhibits a 
jump into a non-uniform one, by opening a "gap" on top of the usual one. As one increases 
Cm the gap continues to grow, as shown in Figures 6,7, and 8 for 5C/(4), SU{5), and SU{10), 
respectively. In particular, we did not observe, in the semi-classical regime, an instability 
towards partial center-symmetry breaking phases, similar to those found in deformed Yang- 
Mills theory and massive QCD(adj), [7,26], see Ref. [27] for a review. 

There is one more interesting issue that appears in the large-A'^c hmit. The density 
of states of large- iVc gauge theories is expected to exhibit an exponential growth, p{E) ~ 
g/3*£; _ ^e/Th ^ where /?* = Th is the Hagedorn temperature. This idea is related to the 
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conjecture that large- A^'c gauge theory is dual to a weakly coupled string theory whose density 
of states is known to grow exponentially. Then, Z{I3) = dE diverges for /3 < /3*, 

indicating the existence of a limiting temperature for the confined phase. As discussed above 
Th is expected to be larger than Tc. In simulations of SU{12) Yang-Mills it was found that 
Tfj = l.llGTc [20]. This seems to imply that it is impossible to extend the confined phase to 
the regime of small-L and weak coupling. Naively, even if one evades deconfinement, it seems 
unlikely that one can evade the Hagedorn instability. It is noteworthy to understand how our 
construction avoids this obstacle. 

The bosonic and fermionic density of states Pb(-E') and pr{E) are properties of the 
Hamiltonian, independent of the boundary condition along the Si. The crucial point is 
that the density operator e~^^ of the thermal compactification is replaced by a Z2-graded 
density operator, e~^^(— 1)^, whose matrix elements, unlike the ordinary density matrix, 
are not necessarily positive definite. The twisted partition function is of the form Z[L) = 
J°°dE {pb{E) - p^{E))e-^^, so that even if the density of states grows exponentially the 
spatially compactified theory need not have a limiting scale. This is most transparent in 
the supersymmetric limit, where all positive energy states are paired bosonic and fermionic 
states, and Z{L) is the supersymmetric index. ^ In the softly broken supersymmetric theory, 
we find a phase transition at the scale given in Eq. (2.22), L* ~ A~^y^ ^ 'XW' ('^^call 
(2.4)), parametrically smaller than the Hagedorn scale f3* ~ A"-*^. Strictly speaking, this is a 
quantum phase transition, which should not be interpreted thermally for small m. However, 
it is continuously connected to the thermal deconfinement transition, and this makes the 
semi-classical analysis feasible. 

3. ^-dependence of the phase transition and topological interference 

In the Minkowski-space formulation of the path integral, the 6-term generates an extra phase 
of the path amplitude, e*'^'^'^) — )• e*'^'''^''^*^'^^'''^ where 7 is the path connecting the initial and 
final field configurations. This is analogous to the Aharonov-Bohm (AB) phase in quantum 
mechanics, and the 0-angle can be viewed as an AB-fiux. Turning on the 9 angle gener- 
ates interference among path amplitudes which affects physical observables and the vacuum 
structure of the theory. 

In the case of Yang Mills theory analytic continuation of the path integral to Euclidean 
space induces a positive definite action for ^ = 0. However, for 6^0, the theory has a sign 
problem, preventing straightforward numerical simulations, see Ref. [28, 29] and references 
therein. In the semi-classical regime the sign problem manifests itself as a complex fugacity 
for topological configurations. For example, the 4d instanton fugacity e~^^ is modified to 
^-Si+t0 _ The density of topological configurations remain unaltered, but they acquire a pure 
phase depending on the topological charge. This leads to qualitative and quantitative changes 

^These considerations are not restricted to supersymmetry theories, or theories with softly broken super- 
symmetry. For example, Nc = 00 non-supersymmetric QCD with multiple adjoint fermions does not possess 
any phase transition upon circle compactification, and satisfies large- TVc volume independence, see [7, 16] and 
references therein. 
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in the vacuum structure of the theory, which we refer to as topological interference effects. In 
particular, the fugacity of 1-defects (monopole-instantons) becomes complex, but the fugacity 
of topologically neutral bions remains unchanged. In the case of magnetic bions the result 
depends on the gauge group G — for SU (Nc) the fugacity remains unaltered, but for G2 it 
acquires a phase depending on the topological charge of the magnetic bion. 

There is a subtle point here. The action of the monopole-instanton is modified as 
exp(— 5*0) — )• exp(— 5*0 + j^) with So = Stt'^ / {g"^ Nc) ■ Since is a periodic variable with 
period 2tt, this implies that there exists a family of Lagrangians: 

= Ckin. + V^^\ k = 0,...,Nc-l (3.1) 

where the Nc potentials Vnp^ are given in Eq. (2.7), and the kinetic term Ckin. in Eq. (2.10). 
For a given value of 6, we must determine the correct effective long-distance theory. The 
correct theory is determined by comparing the vacuum energy density associated with each 
branch, and selecting the branch with the minimum vacuum energy. The theta dependence 
of the vacuum energy density, at leading order in semi-classical expansion, is given by 



£{9) = mm£k{9) = min 



(3.2) 



which arises from monopole-instantons. The contribution due to neutral and magnetic bions 
cancel each other exactly, and the analytic contribution of 4d instantons is exponentially 
suppressed. Eq. (3.2) implies that the vacuum energy is analytic everywhere except for odd 
integer multiples of tt where it exhibits non-analyticity. For example, in the range 6 G (— vr, vr) 
the ground state energy density is £o{9), while for 9 £ (7r,37r) the vacuum energy density is 
£i{9). This implies that at = vr the /c = and k = 1 branches become degenerate. This 
degeneracy that we see in the long-distance effective theory is a manifestation of the CP- 
symmetry in the microscopic theory at = vr. Thus, the vacuum structure and the symmetry 
realization for general group G is modified at = vr as 

-(°) ^ - { f fir — ea P-^) 

For all gauge groups, due to spontaneous breaking of CP, the vacuum in the confined phase 
is two- fold degenerate at ^ = vr. Note that Zg'^ also gives a symmetry that can be used to 
distinguish the phases of G2, F4, Eg at 9 = tt. At very high T, CP is restored and at low T it 
is broken. 

Using Eq. (3.1), we can investigate the 9 dependence for general observables, in particular 
the dependence of the phase transition and the critical scale on 9. It is not hard to see that 
the theta dependence is most pronounced for gauge groups of small rank. We first consider 
the case of SU{2). We will show that L*{9) assumes its maximum at = vr, where the 
monopole-instantons events interfere destructively, resulting in a larger L'^^ for breaking of 
the center symmetry. In the limit where the rank is large, Tc becomes theta independent. 
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3.1 ^-dependence of for Nc = 2 

In the case Nc = 2 the ^-dependence of can be studied quite exphcitly. We take the SU (2) 
roots as ai = —a2 = and find, from the potential given in Eq. (2.7): 



^bion 



2 cosh ( 2V2b') - 2 cos f 2V2a' 



l + i^V2b' 
47r2 



1 



e-^^'cos 



e^"' cos 



V2a' + - + Trk 
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Trk 



(3.4) 



where the various factors of are due to the periodicities given in Eq. (2.2) and the nor- 
mahzation of the roots. Note that k = 0,1 labels the choice of supersymmetric vacuum, 
corresponding to Z4 — )• Z2 i?-symmetry breaking. There is no i?-symmetry when Cm is non- 
zero, and the vacuum branch is the k = one when 9 S (— 7r,7r). For the angular range 
9 G (7r,37r), the vacuum branch is /c = 1. Thus, proceeding with k = and considering 
subcritical Cm, we find that for 9 < tt the global minimum of Eq. (3.4) is at (cr) = 0, {b') = 0. 
For > TT it is located at (a) = {b') = 0, while for 9 = tt the two minima are degenerate, 
and the theory spontaneously violates CP. 

Consider the range < 9 < tt (similar arguments apply to other values of 9). Expanding 
the potential to quadratic order around the minimum (cr) = 0, {b') = gives 
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One of the two eigenvalues of the matrix inside the square brackets in Eq. (3.5) can change 
sign: 
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(3.6) 



where the signs are correlated. Thus, a negative eigenvalue appears at the critical value c* 
Extending the angular range to arbitrary 9, 
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+ ... 



(3.7) 



where ellipsis represent terms with additional (/^-suppression. Note that within the calculable 
regime, the critical Cm given in Eq. (3.7) only acquires 9 dependence due to the ^^.guppressed 
terms. This is special to the = 2 case, for Nc > 3, 9 dependence appear at the leading 
order. 

Using (3.7), we obtain, like vacuum energy density or any other non-perturbative observ- 
able, a multi-branched function, which is two-branched for SU{2), given by 



L*(9) = min 11(9) = L*(0) min 
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(3.8) 
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Both Lq{9) and LKO) are 47r periodic, whereas the min^^o,! ^*k{^) is 27r periodic as it must. 

The physical interpretation of the two branches is as follows. For a given value of 9 the 
potential in Eq. (3.4) always has two extrema. Both of these extrema are center-symmetric 
minima for small g and are located at {a) = 0, For 7^ vr, one of these minima is the true 
vacuum, say /c = and the other is meta-stable, say k = 1. Lq{6) and L^O) are the respective 
transition scales obtained by studying the theory on true vacuum versus meta-stable vacuum. 
At 9 = TT, the two vacua become degenerate, leading to a cusp (non-analyticity) in L*{9), 
where L5(0) = 

As Eq. (3.7) shows, as 9 increases away from zero to vr, the critical decreases. At fixed 
m decreasing implies increasing L^'^ in the range 9 G [0, vr). As in [30], the increase of L'^^ 
with increasing 9 can be attributed to "topological interference". In our context, the physics 
of topological interference is that the effect of monopole-instantons is suppressed at ^ = vr 
with respect to = by destructive interference among path histories. As in Aharonov-Bohm 
phase, this extra topological phase factor gives additional interference among path histories 
(either Euclidean or Minkowski). Therefore, at 9 = ir, the center-destabilizing potential due 
to monopole-instantons is suppressed with respect to 6* = 0. Comparing 9 = n to 9 = it 
requires stronger coupling, meaning a larger value of L, in order for the monopole-instantons 
to overcome the center-stabilizing (^-independent) effect of the neutral bions. 

If the weak coupling behavior holds qualitatively at strong coupling (i.e., for decoupling 
values of the gaugino mass), an increasing L'^^{9) {I3'^^{9)) implies that the deconfinement 
temperature Tc{9) decreases with increasing 9 G [0,7r]. This behavior was discussed in [30] 
and recently observed in lattice simulations [28]. 

4. G2: First order transition without symmetry breaking 

G2 Yang Mills theory is an interesting test case for studying the deconfinement transition 
[31-37]. G2, along with F4 and has no global (center) symmetry and hence no obvious 
order parameter for deconfinement. Lattice calculations show that there is a first order phase 
transition between a low temperature "confined" phase, and a high temperature "deconfined" 
phase. In the confined phase there is an effective string potential at intermediate distances, 
but at very large distance the potential between fundamental charges is screened. The Wilson 
line is small but non-zero in the confined phase, and changes discontinuously at the phase 
transition. In the high temperature phase the Wilson line is close to one.^ 

The group G2 has rank two. The fundamental representation is seven dimensional, and 
the adjoint representation has dimension 14. The tensor product of three 14 dimensional 
representations contains the fundamental representation, [14] x [14] x [14] = [1] -|- [7] -|- . . ., 
which explains why fundamental charges are screened. As in SU{Nc) SYM theory on x 
dynamical abelianization, G2 — ?• C/(l)^, is a consequence of the existence of neutral bions, 

''As discussed above the theory has a global Z2 (CP)symmetry at d = n. This symmetry is broken in the 
confined phase. It is easy to show that it is unbroken in the very high-r deconfined phase. However, it is 
not a priori obvious that CP restoration and the discontinuity of the Wilson line occur at the same critical 
temperature. 
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which generate a repulsion among the eigenvalues of Wilson line, see Fig. 9.^ The long 
distance theory contains the two dual photons ex, their scalar superpartners (j) (the massless 
components of the gauge field along the S^), and fermionic superpartners.^ The fundamental 
Wilson line along is 

Q = e'^-^ . (4.1) 

Similar to the SU (Nc) case, we define the fluctuations b' and a' of the fields cp and a around 
the A;-th supersymmetric ground state 



— log - 

2 167r ^3 



^- ^ (wi log 2-uj2 log 3) + — 6' 



e + 27rk 



27r 



p + 5 



(4.2) 



Here k = 1,2,3,4 corresponds to a choice of vacuum breaking the discrete R symmetry of 
the supersymmetric theory, Zg — t- Z2. One difference from the SU{Nc) case is that the 
expectation value of the fundamental Wilson loop Eq. (4.1) in the supersymmetric ground 
state ((?') = (6) = is nonzero: 



tT{n) = tr e*^'^''^+*fc^'(3'^^°^i"'^i^°s2+a32log3) 

~ tr e '^^-P (^l + i • log ^ - cJi log 2 + 002 log 3^ ^ (4.3) 

^ 5^ /7rlog2 _ 7rlog3 \ _ _q ff^ 
47r V ^ 2 y ■ 47r ■ 



This should not come as a surprise, as there is no center symmetry requiring tr(0) = in the 
confining phase. 

The nonperturbative potential for the b' and a' is due to neutral bions, magnetic bions, 
and monopole-instantons. The result can be derived using the methods introduced in [1], 
based on the non-perturbative superpotential given in [6] and the field redefinition in Eq. (4.2). 



''in phenomenological models it is straightforward to generate eigenvalue repulsion in SU{Nc), but it is not 
easy to find a non-perturbative potential that leads to eigenvalue repulsion in G2, see Section VIII. C in [37]. 
We find that the neutral bion potential automatically gives eigenvalue repulsion for all gauge groups. 

*The dual simple roots al, 02, and the affine root qq = — 2cfi — 02 (two dimensional vectors) have length 
squared 2 (cfi), 6 (02*), and 2 (do), and obey • 02 = —3, al ■ — —1, (S2 • Qo = 0. The fundamental 
weights Uj obey uji ■ a* = Sij {i,j — 1,2). The Weyl vector p = uji + ijj2 satisfies p ■ a*^2 = 1- The 
fundamental representation Cartan generators are normalized as tr HaHt — 2Sab and are explicitly given by 
Hi = 45diag(l, -1, 0, -1, 1, 0, 0), H2 = 4.diag(l, 1, -2, -1, -1, 2, 0). 
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Figure 9: The distribution of the eigenvalues of the Polyakov loop around for a G2 gauge theory, 
for subcritical, critical, and supercritical values of Cm (plotted with ^ = 0.2 in order to visualize the 
jump). 



We find 

hion 

_6e-^'-K+^2) cos{a[ - al) ■ a' - 2e-^'-("S+°i) cos{al - al) ■ a' 
-Cm { (2 + 1^ (2^' • 51 - ^ log 12) ) e-^' "^ cos [a' • a\ + (4.4) 
+ (1 + 1^ • 5^ + J log lOs) ) e-'^'-l cos (^a' . + ^ + ^^^^ 



4 

+ ( 1 + I, (f ■ «1 4 log ^1 1 e-'-"7. cos .35+'' + 

Similar to the SU{Nc) case, we use V^^^^ to denote the overall strength of the bion potential, 
while Cm is the ratio of the strengths of monopole to bion potentials. These parameters are 
given by 

_167rV3 3 6 _jn l_ 

^bion - g2 ^ ^ ^m- ^2^3 ^Igi ' ^ ' 

where ^ ^ 

A3 = e"^ (4.6) 

is the scale parameter of G2 Yang Mills theory, and the explicit factors of g'^ are normalized 
at the scale my\r ~ 1/L. 

Let us now explain the origin of the various terms in Vtotai- Consider first the supersym- 
metric case Cm = 0, when the potential is given by the first two lines in Eq. (4.4). The terms 
independent of the dual photon a' are due to the three kinds of center-stabilizing bions. The 
terms on the second line are due to magnetic bions, which come in two kinds only (because 
03 • a2 = 0). Magnetic bions are responsible for the generation of the mass gap of the dual 
photons and for confinement. It is easy to check that for cj' = 6' = the supersymmetric part 
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Figure 10: Contour plots of the bion- and nionopole-instanton- induced potential Vtotai for G2 for 
subcritical and critical values of Cm (the gaugino mass), as a function of b[ and b'2 (at 9 = 0). Left 
panel: Contour plot of the potential for Cm < < (cm = 1-4, = 1.587). The minimum at 
at the origin ((&') = 0, ((?') = 0) is unique and not destabilized by monopole-instantons. The vacuum 
energy is — 5.6VJ°q„. Right panel: Contour plot of the potential for c^^c^ = 1.587. The minimum at 
the origin is degenerate with the (fo') = (—1.648,-0.345), ((?') = minimum. The vacuum energy is 



of the potential vanishes. It is also straightforward to see that this point is an extremum and 
that the masses of b' and a' excitations are identical, as required by supersymmetry. When 
the gaugino mass is turned on, Cm 7^ 0, the three kinds of monopole-instantons associated 
with the affine and simple roots generate the terms in the potential given on the last three 
lines in Eq. (4.4). These terms can be obtained by a calculation, valid at mL <C 1, virtually 
identical to that of Appendix B of [1], where the SU{2) case was considered. 

Let us now analyze the minima of the total potential in Eq. (4.4). In doing so, we shall 
ignore the 0{g^) terms in the monopole-instanton induced potential. This is consistent in 
the weak-coupling limit. These terms will contribute a shift of the vevs of h' and a' of order 
5^, which can further induce an order g'^ shift in (il), which we ignore. As in the discussion 
of the SU{Nc) case, we first study the expansion of Eq. (4.4) to quadratic order around 
{b') = {a') = and find that a sufficiently large Cm (a number of order unity) destabilizes 
the supersymmetric ground state. It is also easy to see that, as in SU{Nc > 2), there are 
cubic terms in the potential for the b' fields. Thus, one expects a discontinuous transition to a 
ground state where the field b' acquires a vev (and perhaps also a'). To locate the transition 
we study the 9 = case in more detail. In this case the {b') = {a') = vacuum with k = 
has the lowest energy. We find that for Cm < 1.587 the supersymmetric vacuum is 

stable. At Cm = c^, a metastable vacuum with {b[) = —1.648, (62) = —0.345, (<?') = 
becomes degenerate with the zero- vev ground state. Contour plots of the potential are shown 
in Figure 10. 

We find that while there is no symmetry associated with the discontinuous transition to 
a new vacuum at Cm > c^, the value of (tiQ) changes discontinuously from its subcritical 
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value (4.3): 



Cm < cZ : (trf^) = -0.15 ^ , 

47r 

c^ = C+: (trf^) = 3.21 1^ . (4.7) 

To obtain the second number above, we substituted the expectation value of {(p), Eq. (4.2) 
with the vevs of b' as given in the previous paragraph, into Eq. (4.1). It is interesting to 
compare Eq. (4.7) to the results obtained in the lattice simulations of pure Yang Mills G2 
theory [32,34]. Fig. 4 of Pepe et al. [32] shows histograms of (trfi) in the low and high 
temperature phase. The results show that {tril.) changes from slightly negative values below 
Tc to large positive values above Tc, in agreement with Eq. (4.7). 



5. Weak vs. strong coupling non-trivial Wilson line holonomy 

We would like to conclude with a general discussion of the relation between semi-classical 
theories of confinement, discussed in this work, and strong coupling confinement, studied on 
the lattice. For simplicity we consider SU{Nc) gauge theory. A question that is not well 
understood is whether the expectation value of the trace of the Wilson line vanishes in the 
confined phase because 

a ) the Wilson line is dominated by gauge configurations in which its eigenvalues are located 
at the Nc roots of unity with small fluctuations around them. This is the adjoint Higgs 
regime, see Fig. 3b. 

b) fluctuations randomize the eigenvalues over the unit circle, and there is no preferred 
background, as in Fig. 3c. 

This question is a source of confusion especially when one considers the phase transition 
in pure YM theory. There, the transition occurs at the strong scale, hence there is no 
parametric separation of scales to justify an effective field theory in the transition regime.^ 
This regime is often modelled by a potential which breaks the center symmetry in the high 
temperature deconfined phase and restores it in the low temperature confined phase. In the 
limit of asymptotically high T the potential can be justified via a perturbative calculation [3], 
but at low T the coupling is strong and one cannot systematically derive a potential. Ref. [37] 
discusses this issue and proposes that option a) is operative in the low T confined regime of 
Yang-Mills theory. 

First, we emphasize that both a) and b) take place in the confined phase of a large class 
of gauge theories on M'^ x where S]^ is a spatial circle. Examples include A/" = 1 SYM and 

^An exception is the second order transition of pure SU{2) gauge theory. In this case universality arguments 
imply the existence of a 3d effective theory for the Wilson line [38]. We also note that one can always define 
an effective potential for the Wilson line. This potential simply determines the free energy as a function of the 
average Wilson line — it is not the potential in a local effective field theory. 
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QCD(adj), deformed- YM, and deformed-QCD. In all these theories, there are two types of 
confinement operative for different ranges of the parameters. We have 

— - — ^ 1, abelian connnement, ay is operative, 

ZTT 

— ^ 1, non-abelian confinement, b) is operative. (5-1) 

27r 

The small parameter ^^^^ <C 1 arises in the abelian confinement regime as the separation 
of scales between the heaviest dual photon mass vria which is degenerate with given in 
Eq. (2.18) and lightest VF-boson mass m^y: 

rrirr / A \ ^ , A , , 

log «1, (5.2) 



my/ xmw J mw 

implying the first line of Eq. (5.1). This is the regime of abelian confinement where an adjoint 
Higgs mechanism is operative and weak coupling non-trivial holonomy shown in Fig. 3b is 
valid. In this regime, the fluctuations of the eigenvalues are much smaller than the typical 
eigenvalue separation. Once this hierarchy is lost there is no longer an adjoint Higgs effect: 
the eigenvalues have large fluctuations and are essentially randomized over the unit circle. We 
call this the regime of strong coupling non-trivial holonomy, see the discussion in Ref. [19]. 

There is some evidence from lattice simulations that the transition in pure Yang Mills 
theory is between a phase with weak coupling trivial holonomy as in Fig. 3a) and a strong 
coupling non-trivial holonomy shown in Fig. 3c). Fig. 2 of Ref. [39] shows the free energy of 
parameterized Wilson lines O = Diag(e*'^'^, e"*'^'^) in pure gauge SU{2) for different temper- 
atures. At high temperature the free energy agrees with the perturbative one-loop potential 
and has a minimum at = 0. On the other hand, the free energy in the low temperature 
does not have a pronounced minimum at u = ^. Instead, the free energy is essentially u in- 
dependent, suggesting that the Wilson line averages to zero because of the large fluctuations 
of z/ in the confined phase. 

5.1 What changes between the semi-classical and strongly coupled deconfinement 
transitions? 

The phase diagram shown in Fig. 1 contains both types of transitions. The transition from 
the deconfined phase to abelian confinement at small LA continuously evolves to a transition 
to non-abelian confinement at LA ~ 1. Properties of the phase transition that can be studied 
reliably in the abelian regime can be extrapolated to the non-abelian regime. One may wonder 
whether one can improve upon this, and find a more direct connection between deconfinement 
transition in the abelian and non-abelian regime. 

In the semi-classical limit neutral bions Bu = [MiMi] provide repulsion between the 
eigenvalues of the Wilson line. It is not hard to see that magnetic bion-anti-bion config- 
urations, 4-defects, such as [BijBij] = [BijBji] also generate repulsion among eigenvalues, 
and the same is true for 6-defects [BijBjkBki]- In the semi-classical regime these effects are 
exponentially suppressed relative to bions. Recent works provide evidence that these molec- 
ular configurations are semi-classical realizations of IR-renormalons [14,40]. In particular. 
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in the complex Borel-plane, molecular events are associated with singularities located at 
fse.d. ^ ^n, n = (1),2,3, where the existence of the n = 1 term depends on details 
of the particular theory. Note that the location of the singularities is smaller by a factor 
1/Nc relative to the instanton-anti-instanton singularity, and that they survive the large- A'^c 
limit discussed in Section 2.1. The location of the infrared renormalon poles for the theory 
on M^, or for the center-symmetric theory on x E\ in the "^^^ ^ 1 domain, is given 
by tn ~ ^'^i — (1))2,3, . . ., where again, the existence of the n = 1 term depends on 
the theory, and /3o = SA'^c for SYM and /3q = ^iVc for YM. The crucial point is that both 
types of singularities appear with a factor of up to order one pre-factors. The order one 
pre-factors is what seems to be changing as one moves on from the abelian to non-abelian 
confinement regime, and by adiabatic continuity, the IR-renormalons in the strongly coupled 
phase transition replace the role of neutral bions in the semi-classical regime. We believe that 
this direction is worthy of further investigation. 

6. Prospects 

In this work, we provide a microscopic mechanism for the deconfinement phase transition 
in a weakly coupled setting which is continuously connected to the phase transition in pure 
gauge Yang Mills theory with any gauge group. We find that the mechanism for the phase 
transition is universal, independent of the symmetries of the gauge group. Neutral bions 
generate repulsion among eigenvalues of Wilson line, and this effect is counter-acted by the 
perturbative one-loop potential and monopole-instantons which generate an attraction. In 
the weak coupling regime the calculation is based on the semi- classical expansion, and the 
results can be checked (except in the case 0^0) using lattice simulations. It may be possible 
to extend some of the calculations beyond the semi-classical regime using resurgence theory 
and transseries, see [14]. 

The exact counterpart of our analysis can be carried out to describe the deconfinement 
phase transition or rapid crossover in two-dimensional non-linear sigma models, for example 
the 0{N) and CP{N — 1) models, or the chiral principal model. Starting with a supersym- 
metric version of these theories on M x S]^, and turning on a mass term for the fermions, 
a semi-classical transition takes place. We claim that the neutral bions (correlated kink- 
antikink instanton events in that context) are the universal topological configurations leading 
to the confined regime of these theories. 

There are several possible lines of investigation that can be pursued in the future. One is 
the goal to extend our calculations into the strongly coupled regime along the lines discussed 
in the previous section, by studying the singularity structure in the Borel plane both in the 
confined and deconfined phases. Another direction of study is more detailed investigations 
of confinement in the semi-classical regime, for example by constructing effective theories for 
abelian strings, or for the domain walls that separate the different /c-vacua. In addition to 
that we are of course interested in experimental manifestations of the confinement mechanism. 
A recent set of ideas that promises to link deconfinement with the role of topological objects 
is the chiral magnetic effect [41]. This effect leads, in the deconfined phase, to electric charge 
separation in a large magnetic field, in the presence of topological configurations that can 
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generate net chirality. The effect is suppressed by the small density of 4d instantons in 
the high temperature phase, ~ exp[— Stt^/^^]. The density of monopole-instantons scales as 
exp[— 87r^/(5^iVc)], which may result in an enhancement of the chiral magnetic effect. 
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